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O Abstract 

^ The electronic and crystallographic structure of the graphene/Rh(lll) moire lattice is studied 

^ via combination of density-functional theory calculations and scanning tunneling and atomic force 

^ microscopy (STM and AFM). Whereas the principal contrast between hills and valleys observed 

Vh in STM does not depend on the sign of applied bias voltage, the contrast in atomically resolved 

J3 AFM images strongly depends on the frequency shift of the oscillating AFM tip. The obtained 
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I results demonstrate the perspectives of application atomic force microscopy/spectroscopy for the 
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probing of the chemical contrast at the surface. 
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Graphene, a single layer of carbon atoms ordered in a "chicken-wire" lattice [IH3], is 
proposed to be used in many technological applications. Among them are gas sensors [1], 
THz-transistors [5], integrated circuits [6J, touch screens [7], and many others [HE]- One of 
the promising systems on the basis of graphene is its interface with metallic substrates [8]- 
[To] . Here graphene can be used as a protection layer for the underlying substrate [TlHT^ , as 
a spin-filtering material separating two layers of a ferromagnetic material [TSHTTj, or, in case 
of its growth on a lattice mismatched surfaces [for example, Ir(lll), Rh(lll), or Ru(OOOl)], 
as a template for the preparation of ordered arrays of clusters [IHl [IH] • 

In the row of the graphene/metal lattice mismatched systems, the graphene/Rh(lll) 
interface can be considered as an intermediate case between two systematically studied 
graphene-metal systems: graphene/Ir(lll) and graphene/Ru(0001). These three interfaces 
represent the situation of a relatively large lattice mismatch between graphene and a metallic 
substrate. As was demonstrated, graphene on Ir and Ru can be considered as two extreme 
cases of weakly and strongly interacting interfaces, respectively. For graphene on Rh(lll) 
[see Fig. 1(a)] several regions of different arrangements of carbon atoms above a Rh(lll) 
substrate can be found. When using the common notations, the following high-symmetry 
places can be identified: ATOP [A; carbon atoms are placed above Rh(S-l) and Rh(S-2) 
atoms], HOP [H; carbon atoms are placed above Rh(S) and Rh(S-2) atoms], FCC [F; 
carbon atoms are placed above Rh(S) and Rh(S-l) atoms], and BRIDGE [B; Rh(S) atoms 
bridge the carbon atoms]. These places are marked in Fig. 1(a) by circle, down-triangle, 
square, and stars, respectively. Among them, the BRIDGE positions are expected to be the 
most energetically favorable for the nucleation of deposited atoms on top of a graphene layer. 
The nowadays available force spectroscopy and microscopy can shed light on this problem 
and can be used as a tool, which helps to optimize the preparation of ordered arrays of 
clusters on a graphene template [20] . 

In this manuscript we present the combined study of the graphene/Rh(lll) system via 
application of the state-of-the-art density-functional theory (DFT) calculations and scanning 
tunneling and atomic force microscopy (STM and AFM). The calculated imaging contrast for 
STM between all high-symmetry positions for graphene/Rh(lll) is in very good agreement 
with experimental results and this contrast does not depend on the sign of the bias voltage 
applied between a tip and the sample. As opposed to the latter observation, the imaging 
contrast in atomically-resolved AFM measurements depends on the frequency shift of the 



oscillating tip, that can be understood on the basis of measured force-spectroscopy curves. 
The presented results are compared with the available theoretical data. 

The presented in Fig. 1(a) crystallographic model of graphene/Rh(lll) was used in DFT 
calculations, which were carried out using the projector augmented wave method [21], a 
plane wave basis set with a maximum kinetic energy of 400 eV and the PBE-GGA exchange- 
correlation potential [22], as implemented in the VASP program [23]. The long-range van 
der Waals interactions were accounted for by means of a semiempirical DFT-D2 approach 
proposed by Grimme [21]. The studied system is modeled using supercell, which has an 
(11 X 11) lateral periodicity and contains one layer of (12 x 12) graphene on four-layer 
slab of metal atoms. Metallic slab replicas are separated by ca. 18 A in the surface normal 
direction. To avoid interactions between periodic images of the slab, a dipole correction is 
applied [25]. The surface Brillouin zone is sampled with a single fc-point at the F point 
for structure optimization and set to 3 x 3 in the total energy calculations. The STM 
images are calculated using the Tersoff-Hamann formalism [22] , in its most basic formulation, 
approximating the STM tip by an infinitely small point source [23 I2B] • 

The graphene /Rh(lll) system was prepared in ultra-high vacuum station for STM/AFM 
studies according to the recepie described in details in Refs. [19], [29]. The quality and 
homogeneity were verified by means of low-energy electron diffraction (LEED) and STM. 
The STM/AFM images were collected with Aarhus SPM 150 equipped with KolibriSensor 
from SPECS [301 E] with Nanonis Control system. In all measurements the sharp W-tip was 
used which was cleaned in situ via Ar^-sputtering. In presented STM images the tunneling 
bias voltage, Ut, is referenced to the sample and the tunneling current. It, is collected by the 
tip, which is virtually grounded. During the AFM measurements the sensor was oscillating 
with the resonance frequency of /o = 1001541 Hz and the quality factor oi Q = 32323, 
and the frequency shift was used as an input signal in a feedback loop for the topography 
measurements. The oscillation amplitude was set to A = 300 pm. The base vacuum was 
better than 8 x lO^^^mbar during all experiments. All measurements were performed at 
room temperature. 

The DFT-D2 optimized structure is presented in Fig. 1(a) and the variation of the hight 
of the carbon atoms is shown in Fig. SI of the Supplementary material [32]. Carbon atoms in 
the ATOP configuration define a high-lying region sitting at (Iq = 3.15 A above Rh(lll), and 
those in other configurations form a lower region. The buckling in the graphene overlayer 



is 1.07 A. Carbon atoms in the BRIDGE configuration form tlie lowest topograpliic area 
{do = 2.08 A). Tlie HOP and FCC regions are approximately 0.4 A and 0.8 A higher than 
the minima. 

We have estimated influence of dispersion forces on the obtained results: while qualita- 
tively the observed picture remains the same, non inclusion of the van der Waals interactions 
(i.e. standard DFT-PBE treatment) yields larger corrugation {^ 1.8 A) with a very similar 
low region (do = 2. 10 A), but a high region at do = 3.90 A. This is due to the alternat- 
ing "weak" and "strong" interactions of graphene with the Rh(lll) surface. In the case 
of the "strong" interaction between graphene and metal, standard GGA-treatment gives 
reasonable result, whereas for the areas of the "weakly" bonded graphene dispersion forces, 
neglected by the standard procedure, are important. 

Fig. 1(b) shows the DFT-D2-calculated STM image (distance between graphene and the 
tunneling tip is 2; = 2 A) of the graphene layer on Rh(lll). All structural replicas are 
clearly visible in calculated STM image and they are marked by the corresponding symbols 
in Fig. 1(b). The brightest region in this image is the one surrounding the ATOP (A; circle) 
high symmetry arrangement for graphene on Rh(lll). Here the graphene layer is most 
weakly bounded to the Rh substrate. The next two places in the bonding row are FCC 
(F; rectangle) and HCP (H; triangle). The darkest place in the calculated STM image is 
the BRIDCE (B; star) position of the carbon atoms above Rh(lll), where the strongest 
interaction of graphene and Rh(lll) is expected and was calculated. 

Fig. 1(c) and (d) show the experemental STM images of the graphene/Rh(lll) system 
collected over large area and several unit cells of the graphene layer on Rh(lll), respectively. 
The LEED picture of this system is shown as an inset of (c). The results, presented in (d) are 
the part of the results shown in Fig. S2 of the Supplementary material [32]. The presented 
data demonstrate the high quality of the studied system over large regions as well as locally 
and they are in very good agreement with the previously published results [121 [221 E3]- The 
brightest areas are the ATOP (A) positions and the darkest are the BRIDCE (B) ones. 
The clear and unique assignment of all arrangements of carbon atoms in the graphene layer 
and the Rh substrate can be made via comparison of experimental and theoretical data due 
to the extremely good agreement between them (Fig. 1 and Fig. S2 of the Supplementary 
material [32]). The corrugation of the graphene layer on Rh(lll) was measured between 
0.5 A and 1.5 A depending on the imaging conditions. However, the principal contrast in 



STM images is not changed upon variation of the tunnehng conditions (see Fig. S3 in the 
Supplementary material for the corresponding images [S2])- 

The more intriguing results were obtained by AFM imaging of the graphene/Rh(lll) 
system and they are compiled in Fig. 2 and demonstrate the clear atomic contrast in atomic 
force microscopy images of graphene on Rh(lll). In the beginning, the force spectroscopy 
of the studied surface was performed via monitoring the frequency shift of the oscillating tip 
as a function of the distance between the tip and the surface, Af{z). At the same time the 
tunneling current, I{z), was collected. The representative curves for the two most contrast 
areas of the graphene/Rh(lll) topography, namely ATOP and BRIDGE, are shown in 
the upper panel of Fig. 2 as open circles and squares, respectively. The corresponding force 
curves, F{z), are shown in Fig. S4 of the Supplementary material [32]. The tunneling current 
(shown as a semi-logarithmic plot in the inset of the upper panel of Fig. 2) demonstrates the 
exponential dependence on the distance between sample and the tip characteristic for two 
metallic electrodes separated by vacuum. The horizontal dashed lines indicate the frequency 
shifts for the KolibriSensor at which the AFM images were subsequently collected. These 
images with the pointing of the high symmetry stacking positions of graphene/Rh(lll) are 
presented as a low row of Fig. 2 and marked by the corresponding letters. 

The presented force curves can be clearly distinguished from each other that indicates 
the sensitivity of AFM to the different adsorption sites on the graphene layer on Rh(lll). 
This difference in adsorption energies for different places in the graphene/Rh(lll) lattice 
was pointed earlier in Ref. [35]: it was found that the ATOP position is less energetically 
favorable for adsorption and there is a tendency for atoms to nucleate around the BRIDGE 
positions. These theoretical observations are supported by our experimental results: the 
expected equilibrium distance between tip and the sample for the BRIDGE position is 
smaller compared to the one for the ATOP position and the attractive part of the force 
[proportional to the frequency shift, Af{z)] is large for the BRIDGE position. Thus, AFM 
imaging of the surface can be used as a tool for precise detecting of the chemical and the 
electronic contrasts. 

The interesting results are obtained when collecting AFM images at the frequency shifts 
marked in Fig. 2: increasing of the frequency shift during imaging leads to approaching 
the scanning tip closer to the studied surface that yields the increasing imaging contrast. 
This effect is clearly visible in Figs. 2(a-c). In (a) one can see the contrast only between 



"hills" and "valleys" of graphene/Rli(lll), whereas in AFM images shown in (b) and (c) 
the evident difference between ATOP, HOP, FCC, BRIDGE positions (they are marked 
by the corresponding letters) as well as atomically-resolved contrast in the grapehene layer 
are visible. In all these AFM images [Fig. 2(a-c)] the topography of the graphene/Rh(lll) 
lattice as well as all hight variations are well reproduced. For the discussion of the observed 
contrast between HCP, FCC, BRIDGE positions in the obtained AFM images [Fig. 2(d)] 
we can refer to Ref. |34j (see also Figs. S6 and S7 of the Supplementary material [32] for the 
calculated site-projected density of states for carbon atoms and the difference of electron 
density at the graphene/Rh(lll) interface, respectively), where the bonding of atoms of 
different nature on different adsorption sites of the graphene/Rh(lll) is discussed. However, 
the more careful modeling and the further extended discussion is necessary here, that could 
be a topic of the future studies. 

The closer look analysis of Fig. 2(b,c) and Fig. S5 of the Supplementary material [32] 
shows that around "strongly" bonded regions of graphene on Rh(lll) [BRIDGE, HCP, 
FCC) the full unit cell of graphene is imaged when only every second atom in the graphene 
unit cell is visible [larger hexagons in Fig. 2(b,c)]. This is a well known effect for the STM 
or AFM imaging of the graphite- or graphene-based systems. This effect is also clearly 
visible in our STM images (see discussion above and Fig. 1) analogous to the STM imaging 
of the graphene/Ni(lll) system [TTl |35]. Increasing the frequency shift when one goes from 
image (b) to (c) leads to the increasing the imaging contrast and to the appearance of the 
atomically-resolved contrast above ATOP positions, where now every carbon atom in the 
graphene unit cell can be resolved [smaller hexagon in Fig. 2(c)]. This is a first observation 
of the fully resolved carbon ring of the strongly corrugated graphene layer on the metallic 
surface. The present effect can appear due to the local decoupling of graphene from the 
metallic substrate where the local hybridization between electron states of the graphene 
layer and the substrate is very small. A weak hybridization (if any) is also reflected in 
the carbon-atom-projected partial density of states (PDOS) and the difference of electron 
density at the graphene/Rh(lll) interface (see also Figs. S6 and S7 of the Supplementary 
material [32]). The corresponding PDOS for the carbon atoms around the ATOP position 
is very close to the one for the free-standing graphene for which the imaging of every single 
carbon atom in the unit cell was obtained [55^j . 

In conclusion, we present the studies of the graphene/Rh(lll) system by means of DFT, 



STM, AFM, and force/tunneling current spectroscopy. The STM results show the perfect 
agreement with the theoretically calculated images and no principal difference in the contrast 
was found upon changing the sign of the bias voltage between the tip and the sample. On the 
contrary, the atomically-resolved AFM results demonstrate the dependence on the imaging 
conditions, namely on the frequency shift of the sensor. The obtained AFM and STM 
results on graphene/Rh(lll) are understood on the basis of the DFT-D2 calculations and 
they demonstrate the high perspectives of AFM for the imaging of the chemical contrast in 
the heterogeneous systems. 
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Figure captions: 

Fig. 1. (Color online) (a) Crystallographic structure and (b) the corresponding cal- 
culated STM image of graphene/Rh(lll). (c) Large scale STM image of the graphene layer 
on Rh(lll). Tunneling conditions: Ut = +1 V, It = 1 nA. (d) 3D view of the 5 x 5nm^ 
region from (c) showing the structure of the graphene layer on Rh(lll) (see text for detailed 
discussion). Tunneling conditions: Ut = —0.55 V, It = 10 nA. 

Fig. 2. (Color online) (Upper panel) Frequency shift plots as functions of the rela- 
tive distance with respect to the scanning position {Ut = —840 mV, It = 0.48 nA) for 
hills (A) and valleys (B). The inset shows the respective tunneling current, I{z). (Lower 
panel) The AFM images measured at the frequency shifts of the sensor marked by the 
corresponding dashed lines in the upper panel. The large and small hexagons mark the 
graphene unit cell and the carbon ring, respectively. The distorted rhombuses mark the 
graphene moire cell on Rh(lll). 
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Fig. 1, E. N. Voloshina et al. 
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Fig. 2, E. N. Voloshina et al. 
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Supplementary material for manuscript: 
Graphene on Rh(lll): STM vs AFM studies 

E. N. Voloshina/ Yu. S. Dedkov,^ S. Torbriigge,^ A. Thissen,^ and M. Fonin^ 

^Physikalische und Theoretische Cliemie, Freie Universitat Berlin, 14195 Berlin, Ger- 
many 
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List of figures: 

Fig. SI. The height variation of the carbon atoms in the graphene layer on Rh(lll). The 

corresponding high-symmetry positions are marked in the figure. The color bar on the 

right-hand side represents the hight scale in A. 

Fig. S2. Large area atomically resolved experimental STM image (top view and 3D view) 

of graphene/Rh(lll) demonstrating moire structure. The symbols on the left-hand side 

image are the same as in Fig. l(a,b) of the manuscript. 

Fig. S3. Calculated STM images of graphene/Rh(lll) for the distance between the tip 

and the sample of 2; = 2 A. Integration was performed for the valence band stated in the 

energy region oi E — Ep = 0... — 1 eV (left panel) and E — Ep = 0... + 1 eV (right panel). 

Fig. S4. Corresponding force curves, E{z), calculated on the basis of experimental data 

according to the formula presented in J. E. Sader and S. P. Jarvis, Appl. Phys. Lett. 84, 

1801 (2004). 

Fig. S5. AFM image of graphene/Rh(lll): the same as in Fig. 2(c) for better view of the 

effects discussed in the manuscript. Every carbon atom is resolved in the graphene layer 

around the ATOP position. 

Fig. S6. Carbon atom-projected total density of states {a and n) in the valence band for 

the different high-symmetry positions of the graphene/Rh(lll) system. The inset shows 

the corresponding density of states for the Pz character only. 

Fig. S7. Difference electron density, An{r) = %r/Rh(^) — '^Rh(^) — '^gr(^)) plots in units of 

e/A calculated for graphene/Rh(lll). Red (blue) colors indicate regions where the electron 

density increases (decreases). The images are made by means of VESTA visualization 
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software [K. Momma and F. Izumi, "VESTA 3 for three-dimensional visualization of crystal, 
volumetric and morphology data", J. Appl. Crystallogr. 44, 1272 (2011)]. 
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Fig. SI. The height variation of the carbon atoms in the graphene layer on Rh(lll). The 
corresponding high-symmetry positions are marked in the figure. The color bar on the 
right-hand side represents the hight scale in A. 
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Fig. S2. Large area atomically resolved experimental STM image (top view and 3D view) of 
grapliene/Rli(lll) demonstrating moire structure. The symbols on the left-hand side image 
are the same as in Fig. l(a,b) of the manuscript. 
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Fig. S3. Calculated STM images of grapliene/Rli(lll) for the distance between the tip and 
the sample of 2; = 2 A. Integration was performed for the valence band stated in the energy 
region oi E — Ep = 0... — 1 eV (left panel) and E — Ep = 0... + 1 eV (right panel). 
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Fig. S4. Corresponding force curves, F{z), calculated on the basis of experimental data 
according to the formula presented in J. E. Sader and S. P. Jarvis, Appl. Phys. Lett. 84, 
1801 (2004). 
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Fig. S5. AFM image of graphene/Rh(lll): the same as in Fig. 2(c) for better view of the 
effects discussed in the manuscript. Every carbon atom is resolved in the graphene layer 
around the ATOP position. 
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Fig. S6. Carbon atom-projected total density of states {a and vr) in the valence band for 
the different high-symmetry positions of the graphene/Rh(lll) system. The inset shows the 
corresponding density of states for the p^ character only. 
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Fig. S7. Difference electron density, An{r) = ngj./^h{r) — n^^^r) — ng^{r), plots in units of 
e/A calculated for graphene/Rh(lll). Red (blue) colors indicate regions where the elec- 
tron density increases (decreases). The images are made by means of VESTA visualization 
software [K. Momma and F. Izumi, "VESTA 3 for three-dimensional visualization of crystal, 
volumetric and morphology data", J. Appl. Crystallogr. 44, 1272 (2011)]. 
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